Draft Genome of the Leopard Gecko, \u3cem\u3eEublepharis Macularius\u3c/em\u3e

Background Geckos are among the most species-rich reptile groups and the sister clade to all other lizards and snakes. Geckos possess a suite of distinctive characteristics, including adhesive digits, nocturnal activity, hard, calcareous eggshells, and a lack of eyelids. However, one gecko clade, the Eublepharidae, appears to be the exception to most of these ‘rules’ and lacks adhesive toe pads, has eyelids, and lays eggs with soft, leathery eggshells. These differences make eublepharids an important component of any investigation into the underlying genomic innovations contributing to the distinctive phenotypes in ‘typical’ geckos. Findings We report high-depth genome sequencing, assembly, and annotation for a male leopard gecko, Eublepharis macularius (Eublepharidae). Illumina sequence data were generated from seven insert libraries (ranging from 170 to 20 kb), representing a raw sequencing depth of 136X from 303 Gb of data, reduced to 84X and 187 Gb after filtering. The assembled genome of 2.02 Gb was close to the 2.23 Gb estimated by k-mer analysis. Scaffold and contig N50 sizes of 664 and 20 kb, respectively, were compble to the previously published Gekko japonicus genome. Repetitive elements accounted for 42 % of the genome. Gene annotation yielded 24,755 protein-coding genes, of which 93 % were functionally annotated. CEGMA and BUSCO assessment showed that our assembly captured 91 % (225 of 248) of the core eukaryotic genes, and 76 % of vertebrate universal single-copy orthologs. Conclusions Assembly of the leopard gecko genome provides a valuable resource for future comptive genomic studies of geckos and other squamate reptiles

Comparative analysis of evolutionarily conserved motifs of epidermal growth factor receptor 2 (HER2) predicts novel potential therapeutic epitopes

Overexpression of human epidermal growth factor receptor 2 (HER2) is associated with tumor aggressiveness and poor prognosis in breast cancer. With the availability of therapeutic antibodies against HER2, great strides have been made in the clinical management of HER2 overexpressing breast cancer. However, de novo and acquired resistance to these antibodies presents a serious limitation to successful HER2 targeting treatment. The identification of novel epitopes of HER2 that can be used for functional/region-specific blockade could represent a central step in the development of new clinically relevant anti-HER2 antibodies. In the present study, we present a novel computational approach as an auxiliary tool for identification of novel HER2 epitopes. We hypothesized that the structurally and linearly evolutionarily conserved motifs of the extracellular domain of HER2 (ECD HER2) contain potential druggable epitopes/targets. We employed the PROSITE Scan to detect structurally conserved motifs and PRINTS to search for linearly conserved motifs of ECD HER2. We found that the epitopes recognized by trastuzumab and pertuzumab are located in the predicted conserved motifs of ECD HER2, supporting our initial hypothesis. Considering that structurally and linearly conserved motifs can provide functional specific configurations, we propose that by comparing the two types of conserved motifs, additional druggable epitopes/targets in the ECD HER2 protein can be identified, which can be further modified for potential therapeutic application. Thus, this novel computational process for predicting or searching for potential epitopes or key target sites may contribute to epitope-based vaccine and function-selected drug design, especially when x-ray crystal structure protein data is not available

Polymorphism Profile of Nine Short Tandem Repeat Loci in the Han Chinese

Nine short tandem repeat (STR) markers (D3S1358, VWA, FGA, THO1, TPOX, CSFIPO, D5S818, D13S317, and D7S820) and a sex-identification marker (Amelogenin locus) were amplified with multiplex PCR and were genotyped with a four-color fluorescence method in samples from 174 unrelated Han individuals in North China. The allele frequencies, genotype frequencies, heterozygosity, probability of discrimination powers, probability of paternity exclusion and Hardy-Weinberg equilibrium expectations were determined. The results demonstrated that the genotypes at all these STR loci in Han population conform to Hardy-Weinberg equilibrium expectations. The combined discrimination power (DP) was 1.05×10−10 within nine STR loci analyzed and the probability of paternity exclusion (EPP) was 0.9998. The results indicate that these nine STR loci and the Amelogenin locus are useful markers for human identification, paternity and maternity testing and sex determination in forensic sciences

SNP detection for massively parallel whole-genome resequencing

Next-generation massively parallel sequencing technologies provide ultrahigh throughput at two orders of magnitude lower unit cost than capillary Sanger sequencing technology. One of the key applications of next-generation sequencing is studying genetic variation between individuals using whole-genome or target region resequencing. Here, we have developed a consensus-calling and SNP-detection method for sequencing-by-synthesis Illumina Genome Analyzer technology. We designed this method by carefully considering the data quality, alignment, and experimental errors common to this technology. All of this information was integrated into a single quality score for each base under Bayesian theory to measure the accuracy of consensus calling. We tested this methodology using a large-scale human resequencing data set of 363coverage and assembled a high-quality nonrepetitive consensus sequence for 92.25% of the diploid autosomes and 88.07% of the haploid X chromosome. Comparison of the consensus sequence with Illumina human 1M BeadChip genotyped alleles from the same DNA sample showed that 98.6% of the 37,933 genotyped alleles on the X chromosome and 98% of 999,981 genotyped alleles on autosomes were covered at 99.97% and 99.84% consistency, respectively. At a low sequencing depth, we used prior probability of dbSNP alleles and were able to improve coverage of the dbSNP sites significantly as compared to that obtained using a nonimputation model. Our analyses demonstrate that our method has a very low false call rate at any sequencing depth and excellent genome coverage at a high sequencing depth

A survey of overlooked viral infections in biological experiment systems

It is commonly accepted that there are many unknown viruses on the planet. For the known viruses, do we know their prevalence, even in our experimental systems? Here we report a virus survey using recently published small (s)RNA sequencing datasets. The sRNA reads were assembled and contigs were screened for virus homologues against the NCBI nucleotide (nt) database using the BLASTn program. To our surprise, approximately 30% (28 out of 94) of publications had highly scored viral sequences in their datasets. Among them, only two publications reported virus infections. Though viral vectors were used in some of the publications, virus sequences without any identifiable source appeared in more than 20 publications. By determining the distributions of viral reads and the antiviral RNA interference (RNAi) pathways using the sRNA profiles, we showed evidence that many of the viruses identified were indeed infecting and generated host RNAi responses. As virus infections affect many aspects of host molecular biology and metabolism, the presence and impact of viruses needs to be actively investigated in experimental systems

Noninvasive prenatal diagnosis of 21-Hydroxylase deficiency using target capture sequencing of maternal plasma DNA.

Here, we aimed to validate a noninvasive method using capture sequencing for prenatal diagnosis of congenital adrenal hyperplasia due to 21-Hydroxylase deficiency (21-OHD). Noninvasive prenatal diagnosis (NIPD) of 21-OHD was based on 14 plasma samples collected from 12 families, including four plasma sample collected during the first trimester. Targeted capture sequencing was performed using genomic DNA from the parents and child trios to determine the pathogenic and wild-type alleles associated with the haplotypes. Maternal plasma DNA was also sequenced to determine the fetal inheritance of the allele using hidden Markov model-based haplotype linkage analysis. The effect of fetal DNA fraction and sequencing depth on the accuracy of NIPD was investigated. The lower limit of fetal DNA fraction was 2% and the threshold mean sequence depth was 38, suggesting potential advantage if used in early gestation. The CYP21A2 genotype of the fetus was accurately determined in all the 14 plasma samples as early as day 1 and 8 weeks of gestation. Results suggest the accuracy and feasibility of NIPD of 21-OHD using a small target capture region with a low threshold for fetal DNA fraction and sequence depth. Our method is cost-effective and suggests diagnostic applications in clinical practice